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We present the first detection of polarization around the Class low-mass 
protostar LI 157- mm at two different wavelengths. We show polarimetric maps 
at large scales (10" resolution at 350 /im) from the SHARC-II Polarimeter and 
at smaller scales (1.2"-4.5" at 1.3 mm) from the Combined Array for Research 
in Millimeter- wave Astronomy (CARMA). The observations are consistent with 
each other and show inferred magnetic field lines aligned with the outflow. The 
CARMA observations suggest a full hourglass magnetic field morphology centered 
about the core; this is only the second well-defined hourglass detected around a 
low-mass protostar to date. We apply two different methods to CARMA po- 
larimetric observations to estimate the plane-of-sky magnetic field magnitude, 
ON ■ finding values of 1.4 and 3.4 mG. 
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Introduction 



Alignment of du st grains by mag netic fields causes the thermal emission from dust 



to be polarized (e.g., iHildebrandl Il988l ) . The angles of dust polarization measurements at 



millimeter and submillimeter wavelengths are generally considered to be perpendicular to 
the magnetic f ield due to a number of alignment mechanisms (particularly radiative torques, 
Lazarianl 120071 ). allowing for an estimate of the magnetic field direction in the plane of the 

sky. 

Magnetic pressure can support clouds against collapse, though other processes, e.g., 
ambipolar diffusion, can allow gravity to eventually overwhelm magnetic support. The im- 
portance of magne tic fields in sta r formation is not well understood (e.g., weak- versus 
strong- field models, ICrutcherl 120121 ). but observations of the plane-of-sky magnetic field mor- 
phology can provide insight into the coupling of magnetic fields with cores, disks, and out- 
flows. Although flux-freezing during gravitational collapse is expected to create an hourglass 
morphology in the magnetic field lines, the re has only been on e detection of a full hourglass 
around a low-mass protostar to date (e.g., iGirart et al.l 12006). However , seve r al high- mass 
star f ormation regions have clear hourglass detections (e.g., ISchleuningl Il998l ; IGirart et al. 
2009J : iTang et al.lkoioh . 



Large, submillimeter surveys of vario us star-forming clouds and cores have bee n per 



formed with single-dish polarime ters (e.g.. iMatthews et al.ll2009t iDotson et al.ll2010l ). The 



SHARC-II Polarimeter fSHARP. iLi et al.l 120081 ) at the Caltech Submillimeter Observatory, 
is one of the newer submillimeter polarimeters and allows for dual-beam (simultaneous ob- 
servations of horizontal and vertical polarization components) polarimetric measurements. 

Previous millimeter-wavelength interferometric polarimetric maps have been produced 
with the Berkeley-Illinois-Maryland Association (BIMA), t he Owens Valley Radio Obser 
vatory (OVRO), and the Submillimete r Array (SMA) (e.g., IGirart et al.ll2006l ; iKwon et al. 
20061 ; iRao et al.ll2009l ; ITang et al.ll2010l ). but these observations only used single-polarization 
receivers modulated with a quarter-wave plate. The Combined Array for Research in 
Millimeter-wave Astronomy (CARMA) has recently installed 1 mm dual-polarization re- 
ceivers; the first science observations in full-Stokes were made in 2011. 

In this letter we present polarimetric observations from SHARP and CARMA of the 
protostar L1157-mm in the dark cloud L1157. L1157-mm is a low-mass Clas s source (i.e. 
the y oungest of protostars) with a large bipolar outflow spanning ab out 5' (IBachiller et al. 
20011 ) and a perpendicular 2' flattened envelope (jLooney et al.ll2007l ). Though t he distanc e 
to the source is a bit uncertain, with estimates ran ging from 200 to 450 pc (e.g., lKunlll998l ). 
we adopt a distance of 250 pc ( JLooney et al.ll2007l ). 
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These observations trace the magnetic field structure within the core and throughout 
the infall envelope (which, given the adopted distance, could extend to 30-40") of L1157- 
mm. In this letter we will examine how the continuum and the magnetic field morphologies 
correlate with other structures in L1157, such as its outflow and flattened core, and will 
estimate the field strength. 



2. Observations and Data Reduction 

With SHARP and CARMA we create polarimetric maps that show the fractional linear 
polarization, P = \/Q 2 + U 2 /I (where /, Q, and U are Stokes parameters) and the angle 
with respect to the plane of the sky, 9 (measured counterclockwise from north). When 
making uncertainty cuts based on polarization, we used polarized intensity, Pj = IP rather 
than P. Due to the on ly positive nature of P and Pi, de-biasing these values is necessary 



[e.g., 



Vaillancourtl 120061 1 and was done for all observations. 



2.1. SHARP 

Polarimetric observations of L1157-mm were made with SHARP in September of 2008. 
These observations were at 350 /xm (resolution of ~10") and were mosaicked in a chop /no d 



observing mode. The SHARP data reduction is discussed in iChapman et al. I ( jsubmittedl ) . 
The total integration time on L1157 is about 15.8 hours; due to constant chopping, the 
integration time on the source is closer to 7 hours. 



2.2. CARMA 

CARMA observations of L1157 in full-Stokes mode at 1.3 mm were first made in the 
E-array (~4.6" resolution) in July 2011 as a CARMA summer school project. Additional 
observations in the D- (~2") and C- (~0.8") arrays were obtained over th e next year as par t 



of the TADPOL (Telescope Array Doing POLarization) key project (e.g.. lHull et al. 1120131 ). 



An additional D-array summer school observation was also added in June 2012. 

The CARMA dual-polarization system is capable of measuring all 4 cross-correlations 
simultaneously through the circularly polarized feeds a nd full-Stokes corr elator. The reduc- 



tion process for the data in this paper is explained in iHull et al. I ( 120131 ). For L1157-mm, 
the flux calibrator was typically MWC349, but when this source was unavailable, we used 
Mars or Neptune. Flux calibration is accurate within ~15%; however, for the rest of the 
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paper, only statistical uncertainties are discussed. The phase calibrator used for all tracks 
was 1927+739. 



3. Results 

Figure [TJ shows the results (dust continuum and the inferred magnetic field orientation) 
for SHARP and the combined CARMA observations. In the 350 /im SHARP map, an ex- 
tension of the continuum is seen toward the southeast with a positional angle (PA, measured 
counterclockwise from north) of about 124°, differing from the PA of the outflow of 161° 



flBachiller et al. 



in 



Chiang et al. 



2001); however, this extension roughly coincides with N2H + extension seen 



( 12010I ). In the combined CARMA data, the continuum is extended toward 
the northwest with PA m - 20°, which agrees well with ~ -27° from Spitzer 8/xm and with 



-19° from lBachiller et al.l (120011 ) . This extension may be due to photons preferentially escap- 
ing through the poles, thus heating the outflow cavity. The continuum also extends toward 
the east, showing the flattened envelope of L1157-mm; thi s extension matches very well with 
Spitzer and N 2 H + observations of the flattened structure (jChiang et al.ll2010l ). Additionally, 
the continuum observations show an indication of this flattened structure toward the west. 
The east and west continuum extensions are also seen at 1.3 mm with the SMA at ~5.5" 



resolution ( ITobin et al.l 120131 ) . 



For SHARP, all vectors are aligned within ~30° with the outflow angle measured by 
Bachiller et al.l ( 120011 ); however, the southeast vector is offset from the mean of the other 
three vectors by 50°. CARMA vectors have a definitive hourglass shape, with the hourglass 
axis nearly coinciding with that of the outflow (Figure [2]). We note that the difference 
between the hourglass axis and the outflow axis (a proxy for rotation axis) is about 15111 
We do not discuss the reason for this discrepancy, though we note that Kwon et al. (in 
preparation) discusses the possibility of two jets ejected from L1157-mm, and our hourglass 
is tilted toward the more CO-bright (and probably more massive) jet. We note that the 
northwest SHARP vectors also coincides better with the CO-bright jet (within ~15°) than 
the entire outflow (within ~30°). 

As seen in Figure [TJ we also find that for both CARMA and SHARP observations, P 
is significantl y less toward the center of the object, which is typical in polarimetric obser- 
vations (e.g.. iGirart et al.ll2006l ). Several factors may contribute to lower central fractional 



Hull et al. I ( 20131 ) derived a PA difference of 0° for L1157 because t hey defined the outflow direction 



the line connecting the two brightest CO(2-l) peaks, whereas we follow iBachiller et al.l (|200ll) . which uses 
the midline of the outflow cones. 
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Fig. 1. — Polarimetric maps (with polarization vectors rotated by 90° to show the inferred 
magnetic field orientation) of LI 157- mm with the gr ayscale background showing a log-scale 
map of Spitzer 8 /im emission ( iLooney et al.l 120071 ). Magnetic field vectors are shown for 
Pj > la pj and / > 2a /. Left: SHARP continuum contours and vectors are shown in green 
and blue respectively. Cyan shows the 2a intensity detections from CARMA. The length 
of the vectors is proportional to P. SHARP contour levels range from 10% to 90% in 10% 
increments of the peak flux. Right: Cyan (red) contours (vectors) are the combined CARMA 
data at 2.1" resolution and the central blue vector is from SHARP. CARMA contour levels 
are [2, 3,4, 6, 10, 15, 20,40, 60, 100, 140] xa, o=1.02 mJy/bm. Negative contours are not shown. 



polarization, such as: (1) averaging along the line of sight is more likely to smear out the 
polarization through the thickest part of the source, (2) de-polarization at higher density, 
and (3) different grain populations at the center. 



3.1. Comparison between different wavelengths 

To compare SHARP and CARMA values, we smooth the CARMA data to a beam size 
of 10" to provide accurate comparison between wavelengths. We find that the central vectors 
have P=0.7±0.2%, #=-37.9 ±9.0° for SHARP and P=3.80±0.11%, #=-32.2 ±0.8° for the 
CARMA data. These vectors are our only point of comparison, and the angles are consistent 
with each other. This central angle is also consistent at all CARMA resolutions (see Figure 

ED. 



Studies have attempted to explain polarization spectra observed in various molecular 
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Fig. 2. — Hourglass morphology of LI 157 with the red line showing the ax i s of th e hourglass 
and the black line showing the center of the outflow from iBachiller et al.l ( 1200 ll ) . Cyan and 
red vectors are the same as Figure [U Black contours show the CO(2-l) moment-0 map 
(integrated intensity) of the outflow with contour levels of [2,3,4,5,6,7,8,9] xa, a=8 Jy/beam 
km/s (Kwon et al. in preparation). Negative contours are not shown. 



clouds. IVaillancourt fc Matthews! (120121 ) present P[A]/P[350] values at wavelengths, A, of 
60, 100, 850, and 1300 /zm. The median P[A]/P[350] value of a cloud can vary dramatically, 
but the typical value for these wavelengths was ~2. Our value of P[1300]/P[350] « 5.4 
does not fit well within this polarization spectrum. However, comparisons of single dish 
and interferometric fractional polarization should be made with caution; e.g., with varying 
polarization across the source, the interferometer tends to resolve out Stokes / more than 
Stokes Q or U, causing P to be overestimated. Additionally, we only have one point of 
comparison. 
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3.2. Comparison between different size scales 

When a cloud with initially parallel magnetic field lines collapses and flux-freezing 
holds at least partially, the field becomes pinched into an hou rglass morphology. According 



to ambipolar diffusion low-mass sta r formation models (e.g.. iFiedler fc Mouschoviad 11993 
Galli fc Shu! Il993t I Allen et al.l 120031 ) , the hourglass size may range from the size of the pro- 
tostar cloud core (a few hundred to a few thousand AU) to the size of the protostar infall 
envelope (up to about 10,000 AU). Our observations probe resolutions from ~300-2250 AU 
and will provide a diagnostic of the size scale of the pinch. 

From our combined CARMA observations, we apply 3 different Gaussian tapers to 
the visibility data which, along with the non-tapered data, probe 4 different resolutions 
(1.2", 2.1", 3.0", and 4.5") as seen in Figure [3] We note that L1157-mm is only marginally 
resolved at 4.5" resolution. Red vectors are shown when Pj > 2op 7 and are shown spatially 
at approximately Nyquist frequency (~2 vectors per beam). In attempt to see the most 
morphology possible given our sensitivity, we also show yellow vectors for 1.5a p t < Pj < 
2a Pl @ 

Parts of the hourglass morphology can be detected in the lower resolution maps, but the 
full hourglass becomes evident in the two higher resolution plots. SHARP data, however, fail 
to see a definitive pinch even when considering lower a detections (the entire low-a SHARP 
polarization field is not shown in this letter). 

The resolution at which the full hourglass morphology becomes apparent at 1.3 mm is 
approximately 2.1" or 550 AU. We note that this size-scale is highly uncertain (particularly 
due to distance and signal to no ise), though i t is ro ughly consistent with the 290 AU res- 
olution probing the h ourglass in iGirart et al.l ( 120061 ). assuming a distance to NGC 1333 of 
235 pc (ILoinardl l20l3 l . 



3.3. Magnetic Field Strength 

The strength of the magnetic field cannot be measured directly from dust polarization 
maps; however, different methods allow for estimations. The standard technique for cal- 
culating magnetic fields from p olarimetric maps is the Chandraeskar- Fermi (CF) technique 
(IChandrasekhar fc Fermilll953l ) and modifications of the method. 



2 When plotting vectors even down to Pi/crp I =l, we notice that the scatter from the morphology is typi- 
cally small (<10°, indicating we have overestimated errors and/or vectors are not independent). Therefore, 
1.5cr Pj vectors help in displaying the magnetic field morphology. 
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Fig. 3. — The dust continuum emission and magnetic field morphology as measured by 
CARMA to determine the size-scale at which we see the hourglass. Grayscale is the Stokes I 
intensity on a square root scale and the blue ellipses show the beam size for each resolution. 
All vectors have I > 2a j, the red vectors have Pj > 2a p n and the yellow vectors have 
1.5crp 7 < Pj < 1ap I . Top Left: 1.2" resolution, Top Right: 2.1" resolution, Bottom Left: 
3.0" resolution, and Bottom Right: 4.5" resolution, with SHARP vectors for Pi > \.hap I in 
blue (9" resolution). The hourglass pinch becomes the most obvious at higher resolution. 
CARMA contour levels are the same as in Figure 1. The rms error, <r, for each of the four 
panels is 0.94, 1.02, 1.79, and 2.58 mJy/bm respectively. 



We use a modified CF technique (jOstriker et al.ll200ll ) in the large central region where 
polarimetric observations exist (Pi > 2ap 2 and / > 2ai for the combined data at 2.1" reso 
lution). Our application of the technique follows the methodology outlined in iGirart et al. 



( 120061 ). which fits parabolas to the vectors and uses the residuals from the fits to calculate 
the angle dispersion in NGC 1333 IRAS 4A. 

The central region of L1157-mm has a total flux of 0.35 Jy in an area of 70 square arcsec. 
Using a dust temperature of 25 K (Chiang, private communication) and a dust opacity of 
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k„=0.9 cm 2 /g (lOssenkopf fc Henningl Il994l ). we find that the central region has a mean 
volume density of ??,(H 2 )=7.4xl0 6 cm -3 . From N 2 H + observatio ns, we use a non-th ermal 
rms velocity dispersion about the line of sight of 5vi os =0.18 km/s (jChiang et al.ll2010l ). The 
dispersion from the residuals of parabola fits was 7.5° (see Figure H|) and the uncertainty 
of the polarization angle is ~4°, resulting in an angular dispersion of 50=6°. Given these 
parameters, we find a plane-of-sky magnetic field measurement of 1.4 mG. Using these values, 
the mass-to- magnetic flux ratio (using the method discussed in ICrutcherl 119991 ) is found to 
be 1.1 times the critical value for collapse. 
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Fig. 4. — Histogram of angle residuals (difference between measured angle and parabolic fit) 
at each vector location. One vector in the west of LI 157- mm was excluded from the fitting 
since it caused a large reduced y 2 value. 



Since the CF method only calculate s the average m agnetic field strength and has nu- 
merous uncertain ties (see discussion in ICrutcherl 120121 ). we also use another method by 
Koch et al.l ( 120121 ) that calculates the magnetic field strength at all locations where there 
are dust polarization and emission measurements. This technique assumes ideal magneto- 
hydrodynamics (MHD) and that the dust emission gradient indicates the direction of the 
MHD force equation, leading to the equation for the magnetic field strength below: 



B 



sin ip 



sin a 



(VP + pV0) AirR 



(1) 



In this equation, ip is the angle difference between gravity and the intensity gradient, a is 
the difference between polarization and the intensity gradient, P is the hydrostatic pressure, 
is the gravitational potential as a function of radius, and R is the local curvature radius 
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of the polarization vectors. When applying this technique , we assume a densi ty profile of 
the envelope to follow the fits of L1157-mm observations in IChiang et al.l ( 120121 ) (power-law 
of p oc r~ p with p = 2). We also u se a central core mass of M = 0.19 M Q (assumed to be 
similar to L1527. lTobin et al.ll2012f ). Additionally, we assume the pressure gradient, VP, to 
be negligible c o mpar ed to gravity. The local curvature is calculated in the manner suggested 



in 



Koch et all ( 120121 ). 



The magnetic field throughout the large central region where polarimetric observations 
exist (Pi > 2a Pl and I > 2a j for the combined data at 2.1" resolution) can be seen in 
Figure [5j outside this area, magnetic field measurements are unreliable. This figure uses 
median-filtered smoothing at approximately Nyquist sampling (i.e., with pixel sizes of 0.25", 
we replace each pixel value by the median of the surrounding 5x5 box). The colorbar 
indicates the field strength in mG, with an average and median magnetic field throughout the 
region of 3.6 and 3.4 mG respectively. The red contours show the locations where magnetic 
field tension dominates gravity (subcritical). Note that the cen t ral ar eas are supercritical 
(gravity dominates the magnetic field tension). The iKoch et al.l (120121 ) method determines 
the criticality based on sin ip/sin a, while the critical values for the CF method is an average 
about the entire region and depends on calculated values (e.g., magnetic field strength, mass). 
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Fig. 5. — Magnetic field throughout L1157-mm with the color scale in mG. Black contours 
are the same as Figure 1. Red contours indicate subcritical locations, while the rest of the 
colored area is supercritical. 



When using the method from lKoch et al.l (120121 ) to estimate the magnetic field, we have 
3 free parameters: p, M , and distance, d. In the parameter space of our uncertainty in 



- 11 - 



variables, 1.5 < p < 2.5, 0.1 M < M < 0.75 M , and 200 pc < d < 500 pc, the median 
magnetic field follows the approximate relationship: 

px 1 - 9 / M \°'V d ' """' 



When altering these parameters, our median magnetic field value of 3.4 mG does not 
drastically change and will generally be within a factor of ~2. Note that changing these 
parameters does not change the supercritical and subcritical locations in the cloud since 
these locations only depend on the value of sin %p/sm a. 



4. Discussion 

In this letter we present the first detection of magnetic fields in L1157-mm at two 
different wavelengths. Our main results are: 



• 



We find the second instance of a full hourglass morphology in the inferred magnetic 
field around a low-mass protostar. This is also only the second instance around a Class 
source. The axis of the hourglass is nearly aligned with the axis of the outflow, with 
the full detection of the hourglass occurring at a resolution of ~550 AU. 

The angle of the central magnetic field vector agrees at 350 /mi and 1.3 mm and at all 
size-scales. 

We used two methods to calculate the plane-of-sky magnetic field throughout the 
central region and find values of 1.4 and 3.4 mG. 

While this particular hourglass is nearly aligned with the outflow, the TADPOL survey 
results have shown that in general, magnetic field lines are consistent with being preferentiall y 



• 



• 



misaligned (perpendicular) or randomly aligned with respect to outflows (jHull et al. 1120131 ). 



The first well-defined hourglass morphology around a low-mass proto star was observed 



in 1997 in the Class binary system NGC 1333 IRAS 4A (jGirart et al.lll999l ): such a discovery 
has not been published since. Our observation of LI 157- mm is arguably the best example 
of an hourglass shape in any star formation region to date. Characterizing such systems is 
an important step to better understand the role of magnetic fields in star formation. 

IRAS 4 A's hourglass axis of symmetry is misaligned with large scale outflow as measured 



by CO (3-2) (IBlake et al.lll9951 ) by 16° and is over 40° misaligned with the small scale outflow 
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as measured by SiO(l-O) ( IChoil 120051 ). However, SHARP observati ons found that the large 



scale field in IRAS 4A coincides within 1° of the large scale outflow ( lAttard et al.ll2009[ ). Our 
observations of LI 157- mm find that the hourglass axis of symmetry and northwest SHARP 
vectors are similarly misaligned by ~ 15-30°, but the hourglass axis is more aligned with 
the outflow than these SHARP vectors. These results suggest the need to compare more 
polarimetric observations at different resolutions, which will be a future TADPOL project. 

Recent high resolution Expanded Very Large Array observations detect only a single 
continuum source at size scales of ~12 AU (Tobin, private communication), suggesting that 
L1157-mm is most likely a single source. Perhaps cores with a single system are more likely 
to yield an hourglass detection and/or an hourglass axis aligned with the outflow. 

We thank Katherine Rosenfeld, Che-Yu Chen, and Aaron Juarez for their help in pro- 
viding CARMA summer school tracks. 
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